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■ The aim of these notes is to provide a brief review of the topic of neutrino astronomy and in 
I particular of neutrinos from core collapse supernovae. They are addressed to a curious reader, 

■ beginning to work in a multidisciplinary area that involves experimental neutrino physics, 
. I astrophysics, nuclear physics and particle physics phenomenology. After an introduction to 
Q ■ the methods and goals of neutrinos astronomy, we focus on core collapse supernovae, as (one 

I of) the most promising astrophysical source of neutrinos. The first part is organized almost 

• as a tale, the last part is a bit more technical. We discuss the impact of flavor oscillations 
[ on the supernova neutrino signal (=the change of perspective due to recent achievements) 

• and consider one specific example of signal in detail. This shows that effects of oscillations 
I are important, but astrophysical uncertainties should be thought as an essential systematics 

^ ■ for a correct interpretation of future experimental data. Three appendices corroborate the 

, text with further details and some basics on flavor oscillations; but no attempt of a complete 

' bibliographical survey is done (in practice, we selected a few references that we believe are 

I , useful for a 'modern' introduction to the subject and suggest the use of public databases for 

' papers [IJ and for experiments [2 to get a more complete information). 

m : 

' Keywords: Neutrinos, core collapse supernovae, flavor oscillations. 

^ : PACS numbers: 14.60.-z, 23.40.Bw, 26.50.+X, 95.85.Ry, 97.60.-s 

2 ■ Contents 

-4— » ■ 

' ll Neutrino astronomv. methods and goald 1 

' ll.l Main neutrino features! 1 

' ll.2 Concepts of neutrino telescopesi 1 

11.3 Chances for neutrino astronomvl 2 

11.4 Galactic, extraealactic and relic supernovad 3 

12 Supernova neutrino^ 4 

12.1 Gravitational collapse and the 'delayed scenario 1 5 

12.2 Neutrino fluxed 6 

12.3 Effects of neutrino oscillation^ 10 

12.4 Importance of electron neutrino signal (v,, absorption on Argonll 12 

13 Summary and discussioni 14 
IA An example of cross-section (inverse 3 decav^l 15 
IB Fermi integrals and polvlo3 16 
IC A reminder on neutrino masses and oscillations I 17 



1 Neutrino astronomy, methods and goals 



1.1 Main neutrino features 

Neutrinos (and anti-ncutrinos) of electron-, muon- and tau-flavor, are stable, neutral particles. This 
makes them important astrophysical probes; they are expected to point in the direction of the astro- 
physical site of production, as in the more standard case of astronomy with photons.^ Here we have 
in mind the case of 'point astrophysical sources'; but of course 'diffuse sources' are also of importance. 

In normal conditions, neutrinos are invisible. However, they can sometime interact and carry 
away or deposit energy in terrestrial detectors. By contrast, photons are much more easily absorbed 
than neutrinos; they can be observed more easily, but for the same reason their propagation can be 
more easily affected. In certain cases, neutrinos will be the most important signal (think for instance 
to neutrinos from big-bang nucleosynthesis, from the sun, or from a core collapse supernova). 

Some neutrino interactions are of special interest for the following discussion. First, 

z^e + e — + i^e + e, -f e f ^ + e [CC and NC elastic scattering] (1) 

In these reactions an e at rest - say, from an atom - is hit by the neutrino and acquires kinetic energy. 
An important feature is that the hit e maintains the direction of the neutrino when the v energy 
En 3> me ("directionality"). The cross-section is low, dy ~ G\meEv iGp is the Fermi coupling). 
The (lowest energy) neutrino reactions are those of absorption on nuclcons and on nuclei: 

-f J3 ^ e+ + n, i^e + n ^ e + p , , 

z>e + (^,^)^e+ + (A,Z-l), iye + (^,^)^e+(A,Z+l) ^' 

these reactions have usually a threshold, and are only slightly directional (more quantitative statements 
requires care to details, see e.g., App.^. NC cross-sections on nuclei can be as large as ~ Gj^A^E'^, 
when a nucleus composed by A nucleons reacts as a whole (coherent scattering). At higher energies, 
the absorption cross sections on nuclei become ~ Gj^A nipEi, (incoherent scattering). In this case, 
the nucleus is broken and/or hadronic resonances are excited. 

1.2 Concepts of neutrino telescopes 

Let us describe some concepts of neutrino detector, to illustrate what people mean by a 'neutrino 
telescope'.^ (Supernova neutrino detectors fall in the first concept, normally.) 

★ One can instrument a large volume, possibly vetoing for external particle and wait for a charged 
particle coming apparently from nowhere — in actuality, created by a neutrino interaction. (Better to 
be underground for low counting rates, like those related to natural neutrino radiation.) Active volume 
can be a scintillator, a Cerenkov radiator, a layered target, a 'bubble chamber-like' detector. This 
method works from sub-MeV to several GeV energies, because it is subject to the condition that the 
(main part of the) event is contained in the detector. The number of events scales as 

N — Number of targets x cr^ x <I>,y x Time ($ denotes generically a flux). 

In particular, the event rate scales as the volume of the detector. 

^Protons and nuclei of cosmic ray radiation, instead, are deflected by galactic of ~ few and extragalactic 
magnetic fields, at or below nG. They are not expected to point to their sources except perhaps at the very 
highest energies. Fast galactic neutrons instead are another interesting neutral probe. 

^Warning: As it is common in physics, different concepts are blurred and useful at best for orientation; in 
present case, they depend on the type of particle, on the size of the detector ... 
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* One can set a muon counter and timing system underground (or underwater or under-ice) , for muons 
that originate from neutrinos - as those coming from below. Detectors are located underground to 
avoid cosmic ray muons. This is the oldest method and works since muons suffer of mild energy losses 
till 500 GeV (that corresponds roughly to Range^ 1 km in water) . It applies from energies around 
a GeV till several hundred TeV; then the earth becomes opaque even to neutrinos (see e.g., The 
number of events and the i^-induced muon flux scale respectively as: 

N — Surface x $^ x Time where $^ — x a^^ x Na x Range^. 

In particular, the signal scales as the area of the detector (actual target being the earth, the water or 
the ice where the detector is located). 

★ By an extension of previous concept, one could use the earth atmosphere as a target for high energy 
neutrinos to produce inclined air showers; or, use mountains to convert almost horizontal of very 
high energy into visible tau's. In this way, we could observe neutrinos of highest energies. The search 
of inclined air showers is just a spin-off of extensive air shower arrays research activity. Till now, 
however, no positive detection has been claimed. 

In principle we would like to measure a lot of quantities: a) direction of the charged lepton; b) its 
energy; c) its charge; d) tag the flavor; e) tag the time of arrival; f) check occurrence of secondaries 
(n, 7, charged hadrons). In practice, one has to find a compromise between the various and contrasting 
needs of an experiment, e.g. between the wish to have a very 'granular' detector able to see all the 
details of the reaction and the need to monitor a big amount of matter. 

1.3 Chances for neutrino astronomy 

In short, the goal is to use neutrinos to probe astrophysical sources; the information from v can be 
complementary to the one from 7. Some important possibilities in this connection are: 

(1) Solar neutrinos [0.1-20 MeV] There is little doubt that this is V- astronomy'. Among the results 
of a very successful program of observations pioneered by Homestake we quote:'^ a) low energy v 
experiments Gallex/GNO and SAGE prove that the pp-chain (initiated by pp — > De'^Vf.) is the main 
energy source; b) the physics of the center of the sun (pc ^ 150 g/cm'^) is probed. There is consis- 
tency with the theory of solar oscillation eigen- modes (helioseismology) . c) Neutrino oscillations of 
a type predicted in MSW theory are indicated.^ Future observations will aim at the Beryllium line 
(Borexino, KamLAND) and at real time pp-neutrino detection. 

(2) Atmospheric neutrinos [0.05-1000 GeV] primary cosmic rays (CR) come isotropically on earth 
atmosphere and they are not completely understood; they are not thought as astronomy, but they 
belong to astrophysics as much as to particle physics. Atmospheric neutrinos give a very significant 
indication of oscillations, especially thanks to Super-Kamiokande results.^ The study of CR secon- 
daries as the electromagnetic component, muons or atmospheric neutrinos, permits us to investigate 
CR spectra and their interactions with earth atmosphere (which is not that different from possible 
sites of production of CR). In the present context, atmospheric neutrinos will be thought just as an 
important background. 

■^We should recall the important role of certain theorists: J.N. Bahcall, whose activity has been very 
supportive to Homestake since the beginning and G.T. Zatsepin and V.A. Kuzmin, who strongly advocated 
the importance of solar neutrino astronomy. 

^This reconciles SNO observations (1/3 of expected Ve) with those of low energy v experiments (where the 
deficit is less than 1/2). KamLAND experiment supports strongly this picture; more discussion later. 

^MACRO, Soudan2 and K2K support these results. Again K2K, Minos and CNGS long-baseline experi- 
ments will further test these results with man-made neutrino beams. 
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(3) Neutrinos from cosmic sources [unknown energies] This is a vast field and includes a large variety 
of approaches of observation and of objects; presumably, also unknown objects For instance, 
one can search for an excess of neutrino events over the expected background by selecting a solid 
angle-observation window-around a cosmic source (say, an active galactic nucleus) or an appropriate 
time window around a cosmic events (say, a gamma ray burst). Other possibilities are to search for 
self-trigger (excess of multiple 'neutrino' events), or coincidence with other neutrino- or with gravita- 
tional wave-detectors. The observation of point (or diffuse) sources is a very important goal: e.g., v 
(and 7) astronomy above TeV can shed light on the problem of the origin of CR. Till know, several 
experiments hkc LSD, MACRO, LVD, Super-Kamiokande, Soudan2, Baksan, AMANDA, EAS-TOP, 
HIRES and other ones produced upper limits on the fluxes. In future, this type of search will be 
conducted by ANTARES, AUGER, ICECUBE. One of the main hopes is that the neutrino energy 
spectrum remains very hard till ^100 TeV, as suggested by observed gamma spectra at 1-10 TeV 
(another one is that the prompt neutrino background-from charm-is not overwhelming.) 

(4) Supernova neutrinos [few-100 MeV] (this 'cosmic' source is singled out, since it is the topic of the 
rest of the paper). As recalled in next section, most of core collapse supernova energy is carried off 
by neutrinos of all flavors. About 20 events were detected in 1987 by simultaneous observations® of 
Kamiokande II, 1MB, Baksan detectors ^5^ from such a supernova, SN1987A, located in the Large 
Magellanic Cloud, at a distance _D w 50 kpc. Usually, all these events are attributed to inverse f3- 
decay, the one with the largest cross-section (see App.lXll. The experimental detection of these events 
begun extragalactic neutrino astronomy. The agreement with the expectations is reasonable. 

Many operating neutrino detectors like Super-Kamiokande, SNO, LVD, KamLAND, Baksan, AMANDA 
could be blessed by the next galactic supernova. Other detectors like ICARUS and Borexino will also 
be able to contribute to galactic supernovae monitoring in the future. This activity will have a big 
payoff in astro/physics currency: core collapse SN are a source of infrared, visible, X, and 7 radiation 
and possibly of gravitational waves; they are of key importance for origin of galactic CR, for repro- 
cessing of elements, presumably for the dynamics of magnetic fields; they are likely to be related to 
cosmic phenomena like gamma ray bursts; etc. In the following, we focus only on supernova neutrinos. 

1.4 Galactic, extragalactic and relic supernovae 

We close this introduction by classifying and discussing the possible observations of SN neutrinos. 
(Note that, unless said otherwise, the term supernova means always core collapse supernova in these 
notes, even though this is an abuse of notation - supernovae of type la are very important in cosmology 
and astrophysics, and are not core collapse events). 

The hope of existing neutrino telescopes is the explosion of a galactic supernova, for the simple 
fact that the scaling of the flux is severe. In water or scintillator detectors one expects roughly 

300 i>e-events/kton, for a distance D=10 kpc - when our galaxy has a radius of some 15 kpc and we are 
located at 8.5 kpc from its center.^ Various authors estimated the rate of occurrence of core collapse 

®Five other events have been detected by LSD experiment about five hours before the main signal, see 
V.L. Dadykin et at., JETP Lett. 45 (1987) 593. Recently, it was remarked that they could be explained postu- 
lating a pre-collapse phase of emission where only non-thermal of ~ 40 MeV are emitted: see V.S. Imshen- 
nik and O.G. Ryazhskaya, "Rotating coUapsar and a possible interpretation of the LSD neutrino signal from 
SN1987A", to appear in print (preliminary reports presented at 'Markov Readings', INR, May 2003, Moscow 
and LNGS Seminar Series, Sept. 2003, L'Aquila). In this hypothetical phase of emission called also 'cold col- 
lapse' the 200 tons of iron surrounding the LSD detector were the most effective target of terrestrial detectors. 

^One could expect that the chances of getting a supernova where matter is more abundant are higher (the 
galactic center), but one can also object that younger matter, conducive to SN formation, lies elsewhere (in 
the spiral arms). However, we are unaware of the existence of a 'catalog of explosive stars of our galaxy', or 
of calculations of weighted matter distributions of our galaxy. 
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supernovae; for our galaxy, this ranges from ^ 1/(10 y) to ~ 1/(100 y). A recent study of the 
correlations of ^ 200 observed supernovae at cosmological distances with the blue luminosity of their 
host galaxy yields 1/(50-100 y).* A ^ 1/(10 y) lower limit can be already established, since existing 
i/-telescopes did not observe any event yet. Often, one recalls the possibility that SN events can take 
place in optically obscured regions of our galaxy; however, one should also remind that, beside v^s, 
there are other manners to investigate the occurrence of such a phenomenon, e.g., from the released 
infrared radiation. 

Curiously enough, galactic neutrino astronomy is still to begin, but as recalled extragalactic 
neutrino astronomy begun several years ago with SN1987A. In principle, one should profit of the 
wealth of galaxies around us (say, those in the 'local group') to get events at human-scale pace. In 
practice this is difficult, because core collapse SN takes place only in spiral or irregular galaxies and 
not in elliptical ones.^ The only other large spiral galaxy of the local group is Andromeda (M31) but 
(1) its mass is presumably half of our galaxy, (2) its distance is about 700 kpc. A half-a-megaton 
detector (as the one suggested as a followup of Super-Kamiokande to continue proton decay search) 
should get 30 events if efficiency is unit. Perhaps, the best chance would be another SN from Large 
Magellanic cloud (an irregular galaxy) but the odds for such an event are not high. 

Another interesting possibility is the search for relic supernovae, namely the neutrino radiation 
emitted from past supernovae. The practical method is to select an energy window around 20—40 MeV, 
where atmospheric or other neutrino background is small, searching for an accumulation of neutrino 
events there with more-or-less known distribution. The best limit has been obtained by the Super- 
Kamiokande water-Cerenkov experiment |7j, and the sensitivity is approaching the one requested to 
probe interesting theoretical models. In principle, one can suppress the main background (muons 
produced below the Cerenkov threshold) by identifying the neutron from neutrino inverse /3-decay 
reaction. This could be perhaps possible by loading the water with an appropriate nucleus with high 
n-capture rate, that should absorb the neutron and yield visible 7 eventually see e.g., |Hl-^° 



2 Supernova neutrinos 

In Sec l2.1l and 12.21 we present theoretical expectations on supernova neutrinos. More precisely, we 
describe the expected sequence of events of the 'delayed scenario'. This is the current theoretical 
framework '.) 10,, possibly leading to SN explosion. In Sec r2.3l we discuss generalities of SN neutrino 
oscillations. We provide the basic concepts and formulae and discuss the impact on the fluxes. (The 
basic terminology and results are recalled in App. [nj but a real beginner could conveniently consult 
review articles or texts before reading this section. For a more advanced reading, we list in Ref. JM] 
some recent research works on oscillations of supernova neutrinos.) Finally, we complete the discussion 
and show an application of the formalism in Sec l2.4l by considering the reaction Ar— !-K*e^ as a 
signal of supernova neutrinos in an Argon based detector. 

However, the reader should be warned: at present it turns out to be difficult (perhaps impossible) 
to simulate a SN explosion. This could be due to a very complex dynamics; or, it could indicate that 
some ingredient is missing (such as an essential role of rotation, of magnetic fields, etc); or that there 
is nothing like a 'standard explosion'; or, worse, a combination of previous possibilities. In short, we 
have not a 'standard SN model' yet and this makes supernova neutrinos even more interesting. 

*The main unknown comes from the fact that we ignore which is the type of the galaxy that guests us; 
this implies the factor 2 of uncertainty. 

^Their stellar population is older and star forming regions are absent or very rare; in a sense, the stars of 
10-40 solar masses are a problem of youth. 

^"Neutron identification byp-|-n^_D-f7 (2.2 MeV) was proved in scintillators (furthermore, no Cerenkov 
threshold impedes); however no existing scintillator has a mass above 1 kton. 
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2.1 Gravitational collapse and the 'delayed scenario' 



Usually, the life of a star is characterized by a quasi-equilibrium state between gravity and nuclear 
forces. However, the dramatic conclusion the brief-some million years-life of a very massive star of 
~ 10 — 40 Mq is something very different, a core collapse supernova. 

Stellar evolution forms an iron-core, inert to nuclear reactions. This is supported by degeneracy 
pressure of (quasi)free electrons, but when it exceeds the Chandrasekhar mass of ^ 1.4^/© (radius 
~ 3 X 10^ km) it collapses under its weight. The neutron density of the innermost part of the 
core (the 'inner-core', ~ 0.6AfQ) enlarges progressively due to iron photo-dissociation followed by 
electron capture - "infall" phase. When it reaches nuclear densities the increase in matter pressure 
is sufficient to halt the collapse. The 'outer-core' (which is still free-falling onto the center of the 
star) undergoes a bounce on the stiff inner-core. In this moment, an outward-going shock-wave forms, 
producing a prompt neutronization in the shocked material whose mass is about ~ O.4M0 - "flash" 
phase. Then, the shock wave enters a phase of stall, trying to make its way through the outer part of 
the core. This turns the propagating wave into an shock of accretion that involves rest of the initial 
iron core, ^ O.SM© - "accretion" phase. During this phase, convective motions and neutrinos (the 
'delayed mechanism') should revive the shock (that subsequently will eject outer star's layers - the 
SN explosion). The inner core settles in a new quasi-equilibrium state called protoneutron star, that 
smoothly cools and contracts radiating neutrinos of all types - "cooling" phase. Eventually this leads 
to the formation of a neutron star (NS), occasionally seen as a pulsar. Its mass is Mns = 1 — 2 Mq, and 
its radius scales roughly as i?ns ~ 20 km x {Mq/M^sY/^ , due to degenerate character of the equation 
of state. The main features of the collapse process, subdivided in the various phases mentioned above, 
are summarized in Tab.^ 

The most important aspect to note is that the gravitational binding energy released during the 
collapse process (up to the n-star formation) is huge, about 

fs^G^ 7^«(l-5)xl053erg (3) 

(3/5 is for a uniform density distribution) that is about ~ 10% of the n-star rest mass energy M^gC^- 
This is much bigger than the kinetic energy of the ejecta E'kin ^ 10^^ ergs ~ 1% £b (a typical velocity 
of the shock wave is 4-5000 km/s, ejecta mass Moj ~ lOM©). Also much bigger than what is needed to 
dissociate the outer iron core O.6M0/to„ x (2.2 MeV) = 2 x 10^^ erg since the mass of ^^Fe is 123 MeV 
smaller than 13ma -I- 4m„ - but this could be optimistic and the energy losses suffered by the shock 
wave even larger. The energy that goes in photons is very small, i?ium ~ 10*^ erg k. 0.01% £b 
(sufScient to outshine host galaxy though!) and the gravitational wave part is unknown (and depends 
on the detailed dynamics of the collapse) but it is probably even less.^^ The overwhelming part of 
this huge energy is carried away by neutrinos (main reactions leading to v production in the various 
collapse phases are reported in Tab. The neutrino 'luminosity' can be roughly estimated noting 
that ~ 10^'^ erg are emitted in a few seconds in the cooling timcscale, and thus L^, « 3 x 10^^ Lq: the 
supernova neutrino burst outshines the entire visible universe. (Incidentally, we feel there is something 
poetic in these quasi-spherical SN neutrino shells that propagate freely in the Universe). 

^^The gravitational pressure is Pg ~ GjqM^ j . The pressure is Pe ~ ujve, where We = is the 

specific volume and the internal energy u is cpF or pp/{2me) depending on whether electrons are relativistic 
or not (pF=Fermi momentum); thus, Pe ~ he N^^^ or Pe ~ h^/{2me) N^^^ . Since the electron density 
Ne ~ M/{R^m„), non-relativistic e~ lead to the scaling Pe ~ and an equilibrium can be reached; for 

relativistic ones Pe ~ l/R^ and equilibrium is impossible after the core reaches the Chandrasekhar mass of 
M ~ {hc/GNf^^/ml. 

-^^A naive guess is Gn{Mv^ /2)'^ /R ~ EbP"^; it means some billionth of £b with v ~ 4000 km/s. 
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Table 1: Schematic description of collapse and neutrino emission in ttie delayed scenario. (SN progenitor mass 
Mi, ~ 13 ± 3A'/q ). In the first four rows, the main phases are identified. Their conventional names are given in 
Column-1 and the expected dynamics is described in Column-2. Only the main reactions and v-processes ( Column- 
3) are listed. The last row refers to the newborn n-star: when temperature decreases to ~ lO'' K it becomes 
transparent to neutrinos; the emission continues for At ~ 10^ yr, until the temperature drops to ~ 10** K. 
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2.2 Neutrino fluxes 

Here, we describe in some detail the neutrino fluxes. First we discuss the general characteristics and 
present a phenomenological survey, and then we discuss how their luminosity, energy spectrum and 
possible non-thermal effects can be parameterized. We ought to recall the three relevant types of 
neutrino fluxes: 

Ve, Ve and Vx 

where x is anyone among muon and tau (anti)neutrinos. In fact, and have similar properties, 
and Vr are produced by neutral currents (NC) in the same manner and probably, muons are present 
only in the innermost core; thus v^, P^, Vt and Vt should have a very similar distribution. 

Let us begin by describing the general properties of the neutrino fluxes. As seen in Sec l2.2l in the 
delayed scenario the collapse has four main phases. Correspondingly, we distinguish between an early 
neutrino emission, during the "infall" and "flash" phase and a late phase of emission (or 'thermal 
phase'), during the "accretion" and "cooling": see Tab.Hfor more details. The most uncertain phase 
is certainly the one of "accretion", that, together with "cooling", accounts for most of the energetics. 
Perhaps, one could argue that a fair estimate of errors should be just 100 %. In support of this 
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(apparently too conservative) statement, we recall that we have not ab initio calculations of these fluxes 
and alternative (even if incomplete or still speculative) scenarios have been considered. Furthermore, 
the calculations that tried to estimate the effect of rotation (by Imshennik and collaborators and more 
recently by Fryer and Heger [TT] ) found very different fluxes and in particular a severe suppression of 
muon and tau neutrinos. 

Reference ranges on neutrino energies averaged on time (starting at flash time, to = tfi) found 
comparing a number of numerical calculations are: 

(E^^) = 10-12 MeV , 

{EcJ = 11-17 MeV , (4) 
(£;^J = 15-25 MeV , 

The reason of this hierarchy is that neutrinos that interact more - i^e and j>e undergo CC reactions, 
beside NC - decouple in more external regions of the star at lower temperature. In other words, each 
neutrino type has its own 'neutrino-sphere' - v^s one being the outermost. 

The approximate amount (average values over flash and accretion and cooling duration) of the total 
energy £b carried away by the specific flavor is 
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The approximate equality found in numerical calculations has been called 'equipartition', but in our 
understanding, there is no profound reason behind this result. We recall again that these numbers 
should be regarded with caution. 



Next, we would like to introduce a general formalism to describe parameterized neutrino fluxes. 
Such a description requires (a) to know the distance of production D, (h) to assume a distribution over 
the solid angle (usually this is isotropic, up to corrections of the order of ~ 1 — 10 % at most) (c) to 
assign a 'luminosity' function d£i/dt (energy carried by neutrinos per unit time) for each neutrino 
species, (d) to describe the neutrino energy spectrum, presumably black-body of Fermi-Dirac type: 

1 E"^ 

n{E;r]i,Ti) = —- — — — with the normalization factor N = T^F2{r]i) (6) 

iV [1 + exp{E/Tt - ryj)J 

where Ti is the temperature expressed in energy units (the normalization factor N and the meaning 
of F2 are explained in App. Possible non-thermal effects are often described by introducing a 
parameter rji that modifies the shape of the distribution. This parameter is not a chemical potential 
(it is not subject to the condition rj^ = — ryp) and it is often called the pinching factor.^'* (However, it 
is not excluded that the true non-thermal effects are even more dramatic and the high energy tail of 
the spectrum is cutoff as e"*^^/^"-* , where Eq is another new parameter). 
At a distance D from the source, the flux differential in energy and time {t > tfi) is: 

-TErt = TTlir-f^-JT T,) with the normalization factor N' = T,-§Y^ (7) 

dEdt N' AttU'^ dt -t'2[Vi) 



^''if three dimensional effects have an essential role detectable gravitational burst can occur; this adds 
interest in carrying to fulfillment these complex simulations. 

^''This name arises since, for fixed average energy (E), a value 77 > leads to a distribution suppressed at 
low and high energies. The reason 1 10.1 why this happens at high energy is simply that hotter neutrinos are 
in contact with cooler regions than average neutrinos. A typical cross-section that increases fast with energy 
and has a large threshold is Ue^^C ^^Ne~ : changing 77 from zero to 2 decreases by 20 % the event number, 
if (E) = 23 MeV. 
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Figure 1: Fluence spectra forvx-type neutrinos. No oscillation effect is accounted here. The effective temperature 
parameter is set at a reference value T^^ =5.1 MeV. Curves refer to two different values of the effective pinching 
parameter r] — 0,2. The other SN parameters are D — 10 kpc, f — 1/6, £b = 3 ■ IC'^ erg. 



(the index recalls us that we do not take into account oscillations during propagation and the 
normalization factor corresponds to the average neutrino energy N' = (Ei), as recalled in Add. IB)). 
Therefore, one has to calculate or to reconstruct experimentally three functions of the time for each 
type of neutrino, Si, Ti and rji which might be a difficult task. 

For this reason, or just to get a 'synthetic' description, it is common use to introduce the time 
integrated fluxes, i.e. the neutrino 'fluences' from the thermal phase, that are parameterized in a 
very similar manner, namely by (1) an energy fraction parameter fi 

J dt^^ii-Sf,)n£B 

iaccr 

(here, we singled out the energy fraction SfiEs that goes in the 'flash' and fractioned the rest by fi), 
by (2) an effective temperature Ti (time averaged value from taccr) that characterize the spectrum, 
and finally by (3) an effective rji parameter for non-thermal effects. 

In summary, the energy differential fluence, for each neutrino species and for a distance D from the 
source, is given by: 

dF" 1 I E'^ 

1^ = 4^^^- '^■^'^^^^^ jfFU^) [l + exp(£;/T,-^0] 
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Integrating the fluence over the whole surface of emission and over aU neutrino energies 

,.D-^. IdE'^^-^'-'f^^f^'- 



we find the number of i-type neutrinos (N^) emitted during thermal phase (accretion and cooling) 
- oscillations not yet accounted. In Fig. ^ typical ly^-type fluence spectra described by Eq. (jS)) are 
shown. 

We would like to argue that a minimal set of parameters, beside to £b from Eq. J^J and Sfi from 
Tab. n should include the following ones: 

Tp^, K = T^^/Tp^, f and r/ (9) 

These parameters have the following meaning: 

• Tp^ — the effective temperature of electron antineutrinos (presumably, easier to observe); 

UK — increase in temperature of /i/r (anti)i^ (oscillations and NC reactions imply this parameter); 

• f = fu^ = fv^ = the fraction of electron (anti)neutrinos, presumably / = 1/4— 1/6 (see Eq.lSj), 
which constrains f^^ = (1 — 2/)/4 (the case / = 1/6 represents exact 'equipartition'); 

• an effective pinching parameter > 0, equal for all types of neutrinos (that is not expected to 
be accurate, but could be adequate in practice). 

Usually, T^^ is not a very important parameter to describe the neutrino signal, simply because this 
is the lowest temperature, however this can be estimated by a 'reasonable' condition on the emitted 
lepton number — TV? and the parameters of Eq. © : 



T.. = T,J[1 + « - Nl){T,^F,{^)lF,{ij)) I {!£b)\ 

At a further level of refinement, we may introduce time dependent features and distinguish between 
'cooling' and 'accretion' neutrinos. E.g., we have a cooling component whose luminosity dSi/dt scales 
as Tj"*, and whose temperature obeys a time law as: 

T,(t) =r,(0)/(l + t/r) 

(the constant r 10 — 100 sec has to be extracted from the data or computed). On top of that, we 
add for t < Ataccr another rather luminous phase, presumably with a marked non-thermal behavior 
(77 0) and with its own effective temperature. Since the efficiency of energy transfer to matter is 
not large, (anti)i/e should carry a sizable fraction of energy [v^ are of little use to revive the shock, 
but perhaps, only few of them are produced in this phase). 

Sometimes, simplified models of the emission are introduced (see e.g. jOj)- Most commonly, one 
describes the cooling phase as a black-body emission from effective "neutrino radiation" spheres. 
Similarly, one can model the accretion phase by suggesting that the non-thermal neutrino production 
is from interactions with the accreting matter. This suggests that the fluxes are proportional to 
the cross-sections: thus, their scaling should be more similar to E^ than to E^. It is rather interesting 
that there is some hint of such a luminous phase already from SN1987A neutrino signal, see again 
Ref. |12|.^^ In our view, this indication is encouraging for theory and for future observations, even 
though this is not supposed to convince skeptics. 

^^Even if, one could believe that expected deviation from spherical symmetry are large, especially for early 
phases of neutrino emission and for deep layers of the collapsing star. 

^^We would like to comment on the numerical estimate of |I2| . that in SN87A about 20 % of £b was 
emitted during accretion. This is not far from the 'standard' estimate of 10 % reported in Tab.Q however the 
agreement improves further if we assume that i/x are not emitted during accretion (rather than equipartition). 
In fact, A'^^^ (X 20 % £b /6 ^ 0.7 • 10 % fB/2 (the factor 0.7 accounts for oscillations, neglected in Ref. [12\). 
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2.3 Effects of neutrino oscillations 



The basics of three neutrino oscillations and matter enhanced conversion mechanism are briefly set 
out in App. 121 Here we apply them for the supernova, an environment characterized by very high 
electron and baryon densities and of course by very intense neutrino fluxes. 

Let us then start by describing the effect of neutrino oscillations in the stellar medium. Oscillations 
do not affect neutral current events if we postulate to have only 3 types of neutrinos. In fact, the 
fluence + + F° is not changed by reshuffling the fluxes ('NC are flavor blind'). Oscillations 
modify only charge current events (CC). To describe this phenomenon, we need just two functions Pee 
and Pge, the electron neutrino/antineutrino survival probabilities, since the ^ and t flux are supposed 
to be identical. In order to calculate P^e, one has to solve the evolution equation described by the 
effective hamiltonian (see again App. [U)! 

Heff ^ 2.533 ■ t/ ^^^^^'^'^ t/t + 3.868 • 10"'^ pFe • diag(l, 0,0) (10) 
El, 

where H^ff is in m~^, v masses rrii are in eV and the energy E^, is in MeV (similarly for j>e, with 
U ^ U*, and the second term with opposite sign). The first and second term in the r.h.s. of Eq. (|IO|l 
corresponds to the 'vacuum' term and to the 'matter' term, respectively.^* In the latter one the 
supernova density p (in g/cm'^) and the fraction of electrons Ye = Np/ {Np + iV„) must be taken from 
some pre-supernova model. For orientation, a pre-supernova mantle density p = 100 — 200 (ro/r)^ 
g/cm'^ with ro = 10^ km and Yg ^ 1/2 can be used."'^^ 

Inside the core of the star, the 'matter term' dominates and the produced Ve coincides with the 
heaviest state of the effective hamiltonian. It can happen that always coincides with the local mass 
eigenstate during propagation e.g., i/e = ^Ti^) ~ ^^^^ usually called 'adiabatic' conversion. At the 
exit of the star, neutrinos propagate freely as mass eigenstate in vacuum e.g., Ve = ^3, see Fig. |21 [Left]. 
This depends on the unknown size of the vacuum mixing angle 0i3 (C/e3 = sin^ia) and on the electron 
density distribution in the star iVg cx pY^. The approximate values of ^13 when adiabatic conversion 
should occur are shown in the following equation: 

_ I , I ('^\\|2 _ / Ul^^^ if ^13 > 1° adiabatic conversion of j/e —> 1^3 

Pee - Ki^ekeW;i - j ^2^ _ Q 3 if6»^3<o.l° adiabatic convcrsiou of j/g ^ ^ ' 



^'^Indeed, we see a Ue if it stays the same or if or Ur oscillate into u^: Fy^ = PeeF^^ + P^e_F^'^ + PreF^^. 
Rewriting — PeeF^j^ + {Pf^e + PTe)F^_^ and recalling that 1 = Pee + Pfj,e + Pre, we conclude the proof. 

^^There is an additional 'matter term' due to neutrino forward scattering on background neutrinos. Its effect 
on neutrino oscillations with masses as in App.|n|is small, but to see why one needs to go into subtleties. In 
fact, during the most luminous phase (the "flash") the density of background neutrinos A^^ « d£ / dt / {-Kr"^ c- {E)) 
is 10 — 100 X A^e around the point r* ~ 5 • 10* km where ^13 gives MSW conversion. However, the additional 
matter term is strongly suppressed in comparison with usual one, since the relevant current is not Afe(l,0) 
(background electrons at rest) but rather Nvp/E (relativistic background neutrinos). When this current is 
contracted with the current of propagating neutrinos ■u(p')7a(l — 75)ii(p'), it gives zero up to the square of the 
deviation from coUinerity between propagating and background neutrinos, that is ~ (d/r,)^ < 10~^, where d 
is the dimension of the source. The new term is below 1 % at r = r» and scales as r"**; thus its effect is small. 
A strictly related and more detailed discussion is in Y.Z. Qian and G. Fuller, Phys.Rev.D 51 (1995) 1479. 

^^Note that we are assuming that the pre-supernova dynamics does not modify in an essential way the 
structure of the mantle of the star. The modifications due to the shock wave, usually, do not lead to large 
effects. However it is possible at least in principle that a massive occurrence of stellar winds, explosive nuclear 
reactions, and/or instabilities modifies the mantle before the occurrence of the core coUapse. These possibilities 
will be better tested by astronomical observations of the pre-supernova, by the study of the supernova spectra 
and/or possibly by theoretical modelling of the star mantle; e.g., using the delay between the neutrino burst 
and the light. (We thank Marco Selvi for this important remark.) 
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Figure 2: These plots show the effective neutrino masses of neutrinos [Left panel] and antineutrinos [Right 
panel] inside the star. As visible from r ^ oo regions, we assume that the neutrino spectrum obeys a 'normal' 
mass hierarchy. (We do not emphasize another possibility compatible with the data, 'inverse' mass hierarchy. The 
common mass scale is immaterial for oscillation however.) 



(at present, we cannot exclude that 9i3 falls in an intermediate case).^" Solar neutrino mixing makes 
almost certainly adiabatic the second conversion, if the first should fail. 
Thus, accounting fo oscillations, the fluence of Ve becomes: 

F -/ if^i3>l° .,2^ 

" \ O.SF^^ + 0.7F°^ if 013 < 0.1° ^ > 

Similarly, for antineutrinos i>e — i>™{t), see Fig. 12 [Right], that implies Pgg = | (j/e|i'e(i)) P — U^i ^ 0.7. 
Thus, the formula for the flux becomes Fp^ = 0.7Fl^^ + 0.3F°^. Now we can make the argument for 
oscillations: Since we expect that F^^^ ^ and F^^ ^ , oscillations should modify the expected 
supernova neutrinos fluxes. These modifications are large (e.g., the flash yields little in CC: NC events 
range from 70 to 100 %) and can be observable, but the message that we want to stress here is sim- 
ply that these effects should be taken into account in order to interpret the SN neutrino signal correctly. 



^"This qualitative discussion of neutrino oscillations in matter as illustrated in Fig. |21 corresponds to the 
approximated analytical expression for the probability of survival Pee ~ Ue2PH + U^ail — Ph), where the 'flip 
probability' associated with is Ph = exp[-f/f3/(4 • 10"^) • (20 MeV/S^)^''^]. (The analytical expression of 
the exponent is 27rro(\/2GFiVeo)^''^(Am^/2£;)^/^ and assumes that Ne{r) — Neo {ro/r) , see last reference in 
|13|'). The corresponding flip probability associated with solar mixing is Pl = with good approximation. 
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Finally, we consider the 'earth matter effect', possible operative if SN neutrinos cross the earth 
before hitting the detector. We will show that, with the current oscillation parameters, it is not very 
large. As we saw, in a possible scenario (=normal mass hierarchy, very small ^13) neutrinos exit from 
the star as {ve) — > W2) and — > due to the MSW effect [T^, or in other words. 



Pee = sin^ 6*12 ~ 0.3 and Pgg — cos^ 612 ^ 0.7 



If (anti)ncutrinos cross the earth in the last stage of their path, new oscillations will occur (since 
vacuum eigenstates are not eigenstates in the earth matter) and previous expressions will be modified. 
For constant density (say, earth mantle - « 4 g/cm'^) the solution of a two-flavor version of Eq. (|10|l 
gives: 



Pee = sin^ 612 



4ecos2 6'i2 . 2 ( ^'^12^ nz 7^5 ^ ttt- 

1 + /-, , -A TH- ■ V 1 + -4£cos2 6li2 

(1 + e)2 - 4ecos2 6*12 \ AE 



(13) 



with 6*12 « 33° and Ato?2 ~ 7 • 10"^ eV^, where 



^ V2Gf iVee pe/(4 g/cm^) ■ 1^^/(0.5) ■ E/{2Q MeV) 

^ ^mlJ2E - ^" ^mlJ{7 ■ 10-5 eV^) 

For Ve, just replace 612 90° — 6*12. Earth matter effect is larger than for solar neutrinos, simply 
because supernova energies are larger, see Eq. ((117)1 . This can give rise to spectacular wiggles, especially 
if large energies events are seen. Numerical considerations based on previous formulae suggest that 
this investigation will be demanding. If (or when) the position of the supernova will be known, it will 
be possible to include such an effect, reducing ambiguities in the interpretation of the signal. 



2.4 Importance of electron neutrino signal (z/g absorption on Argon) 

In order to complete the discussion and to show an application of the formalism, wc will consider in 
detail the specific supernova neutrino signal provided by the reaction of absorption 

z^e + Ar ^ K* + (14) 

that has a large cross-section. The signature for reaction (|14|l is given by a leading electron accompa- 
nied by soft electrons from conversion of K* de-excitation 7's in the Argon volume surrounding the 
interaction vertex. This signal could be seen by the forthcoming detector ICARUS ^B] based on the 
liquid Argon technology.^^ (We are not going to discuss the more difficult and important question of 
'what we can learn from supernova neutrinos', whose answer will of course depend on which neutrino 
detectors will be working when next galactic supernova will explode and what will be the distance of 
this supernova; but it is almost from granted that we will learn a lot from the VeP — > ne^ reaction for 
the reasons recalled in Add. fX|) . 

^^In this discussion, we want to emphasize the potential of an ideal i/^ detector, putting aside technical 
limitations like need of a long term stability of operation, detection threshold, efficiency and finite resolution. 
Let us recall that also other detectors can see the z/g signal with other reactions, even if usually this is not 
the main signal. For instance: (1) + ^^O ^ e + ^^F (with Q — 15.4 MeV) can be exploited at water 
Cerenkov detectors as Super-Kamiokande or SNO due to the angular distribution (^^F rapidly decays by 
proton emission); (2) Ve + ^^C e 4- ^^N (with Q — VIA MeV) can be seen in scintillators detectors (LVD, 
Borexino, KamLAND, BAKSAN), with the great advantage of offering a double tag, due to the (3'^ decay of 
Nitrogen; (3) + O ^ e + p + p (with Q = 1.4 MeV and a large cross-section) can be used at the inner 
part of SNO (the signal is given by a lone electron, in contrast with neutral current, or electron antineutrino 
reactions on deuterium that are tagged by additional neutrons). detection profits of staying closer to the 
philosophy of solar neutrinos and of employing literally solar neutrino detectors. Note that a big Q value or 
a rapid rise of the cross section amplifies the difference between the case with and without oscillations. 
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In a 3 kton liquid Argon detector the number of i/g-absorption events is about 400, for 

a supernova exploding at D = 10 kpc. To calculate this number, one simply multiplies the fluence 
(including oscillations) by the number of target nuclei and by the cross-section of the reaction, and 
than integrates over the possible neutrino energies. In the present calculation we employed a 'hybrid 
model' for the cross-section of the reaction in Eq. H14() : shell model for allowed transitions, and 
random phase approximation (RPA) for forbidden ones (see below). The other inputs were: (a) a, 
normal hierarchy of neutrino masses; (b) 613 large enough to produce Pee ~ - that is, Ve 1^3 
producing F,y^ = ; (c) a,n exact equipartition of the fluxes (/ = 1/6) and = 3 x 10^^ erg; (d) a, 
spectrum without pinching (77 = 0); (e) sl neutrino temperature of Ti^^ =5.1 MeV, corresponding to 
an average energy = 16 MeV (following the indications of the most recent calculations (5], we 

assume that in absence of oscillations the temperature of electron neutrinos is T^^ =3.5 MeV, that is 
closer to T^^ than thought in the past). 

The use of adequate cross-section for neutrino absorption reaction on Argon is important. Allowed 
transitions to low lying Potassium (K) excited levels^^ dominate for neutrino energies less than ~ 
15 MeV (i.e. in the energy range of interest for solar neutrino experiment). Shell model computation 
[T7| allows to reliably describe the allowed cross-section. At higher energies, as for the SN case here 
considered, forbidden transitions become relevant as well. These are dominated by the collective 
response to giant resonance, so that the RPA model ^] is usually considered sufficient to describe 
the non-allowed contributions to the {i/e, Ar) cross-section.^'^ 

How the number of events changes, with reasonable changes of the input parameters? To answer 
this question, we can calculate the percentage variation 100 x SN/N of the number of absorbed 
under a number of alternative hypotheses: 



T + AT 


T-AT 


/^l/8 


,7^2 


Pee 0.3 


+51 % 


-45 % 


-25 % 


+ 15 % 


-16 % 



The first two columns show the effect of changing the temperature by AT = ±1.3 MeV; the third 
column, describes the effect of non-equipartitioned fluxes; the fourth one, the effect of having a pinched 
('non-thermal') spectrum; the last column, assumes that z^e ^ 1^2 due to very small 613. This shows 
that the present uncertainty in the temperature has a big impact on the expected signal, about 50 %. 
It shows also that a mixture of various phenomena can affect the flux at the ^ 20 % level. To 
separate these effects clearly, it will be important to study several properties of the neutrino signal, 
like distributions in time and energy and use several reactions. In Fig. El we show the calculated 
number of expected events for a wide range of values of the effective temperature. 

In some situations, the electron neutrino signal can lead to 'model independent' inferences. For 
instance, if it were possible to demonstrate that the earth matter effect (associated with solar Am^ « 
7.1 • 10^^ eV^) occurs in ly^ events and in events, we would have a proof that 6*13 is small. If instead 
it occurs only for D^., the converse is true and furthermore, the hierarchy must be normal (there 
is an adiabatic conversion associated with the heaviest neutrino). It should be remarked however 
that a 'golden' observation (that is seeing one or more wiggles) requires a great precision in energy 
measurement or a lucky configuration, namely, a supernova exploding just below the horizon. In fact, 
the phase of oscillation with solar Am^ is close to 7r/2 for lengths of propagation through the earth 
of the order of 350 km x £:^/(20MeV), see Eq. ((T^. 

^^The allowed transitions in Eq. 1141 include two contributions: (1) Fermi transitions from ''"Ar (J'^ — 
0+,T = 2) to the isobaric analog state of ""K (J'' = 0+,r = 2) at an excitation energy of 4.38 MeV and 
(2) Gamow- Teller transitions to several low lying J'^ — 1"^ ,T = 1 states of with excitation energies 
between 2.29 to 4.79 MeV. 

■^''in the RPA calculation of Ref. |18j, all forbidden transitions to '"^K levels with J < 6 and both parities 
have been included. 
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Figure 3: Number of expected Ve events as a function of the effective temperature. Oscillation are separately 
accounted for - full line and dotted line - according to the two reference cases of Eq. <J2I) . Neutrino fluence is taken 
from Eq. with T free to vary and pinching parameter set at 77 = 0. The other SN parameters are: D — 50 kpc, 
/ = 1/6 (strict equipartition), and £b ~ 3x 10^'^ erg. The cross-section 'hybrid model' for Ve-absorption on Argon 
is used. A reference 3 kton detector mass (corresponding to 3 x 1.5 ■ 10"^^ Ar targets) is considered. For simplicity, 
we assume an ideal detector, without threshold on the final state electron energy and full detection efficiency. 



But note that even the absence of a i>e signal would be a precious information. Indeed, to help 
the 'delayed explosion' to take place, it would be better to have a depletion of during accretion. In 
that case, the number of i^e events in the first half-a-second should be small, due to Eq. (|12|l . whereas 
they should be seen during cooling. Similarly, if non-standard scenarios (like collapse with rotation) 
are realized, Vx can be depleted also during the cooling phase. In this case, Ve events would be rare 
even during cooling. 



3 Summary and discussion 



In this introduction to i/-astronomy, we focused mostly on supernova neutrinos. We aimed at helping 
the orientation of a reader in this field, so we did not attempt to give a comprehensive study (i.e. we 
did not consider all theoretical possibilities or scenarios, or reactions to detect neutrinos). Rather, 
we offered a selection of the background information, provided some few formulae, reference numbers, 
and showed illustrative calculations. Let us conclude by recalling some of the important points we 
touched: 
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★ Neutrino astronomy is theoretically appealing and rich of promises. Supernova neutrinos are a very 
well defined and interesting possibility. 

★ Neutrino observations from SN1987A are not in contradiction with the general theoretical picture. 
However, supernova explosions are still mysterious, and this warrants more discussion and stimulates 
more efforts. 

★ Next galactic supernova will permit us much more precise observations and this will be certainly very 
helpful to progress. In particular, the response from new generation neutrino detector(s), sensitive to 
various types of neutrinos and reactions could be of major importance (we discussed in some detail 
the case of a detector like ICARUS, that combines a large mass with a high resolution and detection 
efficiency). 

★ The effects of oscillations are important and have to be included. Conversely, one could combine 
experiments and use theoretical information in order to attempt to make inferences on oscillations, 
but astrophysical uncertainties should be thought as an essential systematics for this purpose. (In 
other words, there are chances to learn something on neutrinos, but, in our view, the primary aim of 
these observations is just supernova astrophysics.) 

★ All this is fine; the most important task left is an exercise of patience, — 100 years for next galactic 
supernova. 
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Appendices 

A An example of cross-section (inverse /3 decay) 

The 'inverse /3-decay' reaction 9e +p — s- e+ + n is particularly important for actual neutrino detection. 
Indeed, it has a large cross-section and water Cerenkov and scintillator based detectors have many 
free protons. For illustration, we recall here a simple approximation of this reaction (from last paper 
of Ref. ^3]) and refer to ^H] for a more complete discussion. 

The tree level cross-section in terms of Mandelstam invariants s, t, u is: 



da Gjp cos^ 6c 
dt 27r(s — TOp) 



[A{t) - B{t){s -u) + C{t){s - u)2] (A.l) 



A, B, C are well approximated as ci 4- 1 • C2 at the energy of supernova neutrino detection: 
A w AP{1 - g^){t - ml) - lVPl^^{l + g^) - 2mlMAg{l + 

B^g{l + Ot (A.2) 
C^{l + g^)/A 

where M = (to„ + mp)/2, A = 77i„ - m^, ^ = 3.706 and g = -1.270 ± 0.003. Eq. (fO|l is related by 
Jacobians to the cross-sections differential in the lepton energy E,, , or in the angle 6 — ~ t: between 
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the incoming neutrino and the charged positron: 

da da da epe da 

—— — zmr,— — and — = p -— (A. 3) 

dE^ ^ dt dcose 1 + e(l _ |i cos 9) dE^ ^ ' 

Of course, for the first formula one has to express the Mandelstam variables in terms of E^ and E^^^ 
e.g., t = rn^ — — 2nip{Eiy — E^)- To evaluate the second formula, one first defines e = E^/nip and 
calculates the positron energy E^ (and momentum pe) from E^ = [{E^, — S){1 + e) + ecosO{{E^ — 
(5)2 - m^fc2)i/2]/^^ where 5 = {ml - - ml)/{2mp) and fc = (1 + e)^ - (ecos6')2. Note that E„ is 
one-to-one with Ee at zeroth order in e. The threshold of the reaction is at E^ > 1.806 MeV. 



B Fermi integrals and poly log 

Let us consider the function fn{x, rj) where n = 1,2, 3... and 77 is a real parameter: 

a;" 

f„ix,r]) = with a: > 

1 + ^ 

This function is needed to define the Fermi integral of n-th order as follows: 

/"OO 

Fn{ri) = / /„(x,77) dx = -Li^+il-e") -n! (B.l) 



The last expression involves the polylogarithm function Li„(a;). 

• This integral appears commonly when using the Fermi-Dirac distribution, that can be written as 
T^f2{x, rj), with x — E/T. E.g., integrating this distribution for all values of the energy E we get the 
normalization factor of Eq. jnj of Sec 12. 21 

/•OC 

/ T^f2{x,v) dE = F2{v) (B.2) 
Jq 



• Eq. (|B.1(I is also useful to express energy momenta: 

^^^^^'W)' ^^'^^^'■fH' (B.3) 

The variance {SE) stays constant at better than 0.2 % for 77 < 5 at the value (SE) ~ 1.73 -T. A useful 
approximate expression for the average energy (SE) can be obtained using ^{^{—l+y) /Li^l—l+y) — ! w 
(50.5 — 41.8?/ — 6.4y2 — 2.3?/'^)/1000, where we have in mind the identification y = 1 — exp(— 77). 

• For ?7 = 0, 1, 2 we get F2 = 1.803, 4.328, 9.513 and F3 = 5.682, 14.39, 34.30. 

• The series expansion Li„(z) = J2m=i -z™/™" of the polylog leads to some identities: 



Lio(z) = z/(l - z); Lii(z) = - log(l - z); Li„+i(z) 



/■^ dr 

/ Li„(T)— with z > (B.l) 
Jo T 



At z = ±1, the polylogarithm can be expressed by the Z-function (See Eq. HB.ll) '): 

Li„(l) ^ Z„, (^2,3,4,5 = 1.645, 1.202, 1.082, 1.037) 

Li„(— 1) = —(1 — 2^~")Z„ (connection of Bose and Fermi integrals) 



(B.5) 
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C A reminder on neutrino masses and oscillations 



(1) The mixing matrix U (introduced by Sakata and collaborators in 1962) connects neutrino fields 
of given flavor and of given mass: 

vi{x) — Uii Vi{x), where I — e, and z = 1, 2, 3 (C.l) 

This implies a relation between states of ultrarelativistic neutrinos. In fact, the decomposition in 
oscillators v{x) ^^^^('^pAWpAe'^^ + b\xVpxe~'''P'^) imphes b\ = Uuh] and a\ = Ul^a\, so: 

No change if the type of mass is Dirac instead than Majorana (last one being theorist's favorite). 

(2) From previous considerations, it follows that if we produce a state of flavor ^ at t = it will 
acquire overlap with other states at later time ( "appearance" of a new flavor) and at the same time 
it will loose overlap with itself ("disappearance"). This was shown by B.Pontecorvo in 1967, though, 
the first idea dates back to 1957. Thus, a state with momentum p becomes 

We{t)) = C/;>.(i)>, where \,,,{t)) = e'(P-E.t) |j.^(o)) (C.3) 

The energy of neutrinos with different masses cannot remain the same in the course of the propagation, 
since Ei Ki p + mf j {2p) (ultrarelativistic approximation always applies to the cases of interest). When 
the distance between production and detection satisfies L ^ t ^ E/Arnfj, \ve{t)) becomes different 
from \i'((0)), if the mixings Uu are large enough. As usual, Amfj = rrij — m|. 

(3) The effective hamiltonian of propagation in vacuum is Udia,g{m^)W /2E (for antineutrinos, U — > 
U*), but in matter there is an additional term ±V2G'FAfediag(l, 0, 0) (+ is for Ve, — for De] Ne = 
pYe/nin — number density). This term is linear in the Fermi coupling. It describes a coherent 
interaction of neutrinos with the matter where it propagates. This can drive neutrinos to be "local" 
mass eigenstates during the propagation, thus exiting from a star as vacuum eigenstates, i.e.: \ve) 
\i'2} or Ive) \i^3)- In the sun, this effect (named after MSW after Wolfenstein, Mikheyev and Smirnov 
|13|) is partial and it is pronounced for highest energy neutrino events, e.g., the CC events at SNO. 
In the supernova, it can be complete. 

(4) Putting aside the indications of LSND indication (that will be tested at MiniBooNE) we know 
from a number of experiments that the usual 3 neutrino flavors most likely oscillate among them and 
this points to the following (roughly 1 sigma) ranges of the parameters |19| : 



= 7.1 ±0.7- 10-5 eV^; 6I12 = 33° ± 2° 
t23 = 2.0 ±0.4- 10^3 gv2. 6(23 = 45° ±7°; U13 



AtoC = 2.0 ±0.4- 10^3 gv2. = 45°±7°; 6*1 3 < 9° ^^''^^ 



The 3 mixing angles given above parameterize the unitary mixing matrix Uif. 

\Ue3\ = sin6'i3, \Ue2/Uei\ = tan 6'i2 (solar mix.), lU^^/Ural = tan6'23 (atmospheric mix.). 

(5) We have some bounds on neutrino masses (in this sense, we know something more than the 
oscillation parameters Am^j) from other sources: from /3-decay (Mainz, Troitsk), v^^TWSl™? — 
2.2 eV; from neutrinoless double beta decay (Heidelberg-Moscow at Gran Sasso, IGEX), | J^i Ulimi\ < 
0.3 — 1 eV [a claim was made that the transition has been observed, but in our opinion, with a weak 
significance]; from galaxy surveys (2dF) ™« < 1.8 eV or combined cosmological results (including 
WMAP) ^ ■ nii < 0.7 eV. We do not discuss them further, but we note that they suggest that the 
kinematic search of effects of neutrino masses with SN neutrinos is difficult or impossible. 
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